Regimes of strong light-matter coupling under incoherent excitation 
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We study a two-level system (atom, superconducting qubit or quantum dot) strongly coupled to 
the single photonic mode of a cavity, in the presence of incoherent pumping and including detuning 
and dephasing. This system displays a striking quantum to classical transition. On the grounds 
of several approximations that reproduce to various degrees exact results obtained numerically, we 
separate five regimes of operations, that we term "linear", "quantum", "lasing", "quenching" and 
"thermal". In the fully quantized picture, the lasing regime arises as a condensation of dressed 
states and manifests itself as a Mollow triplet structure in the direct emitter photoluminescence 
spectrum, which embeds fundamental features of the full-field quantization description of light- 
matter interactions. 



I. INTRODUCTION 

The strong coupling regime is the ultimate limit of 
light-matter interaction, at the level of a single quantum 
or a few quanta of excitations. This gave rise to the field 
of cavity Quantum Electrodynamics (cavity QED) [I], 
which in the recent years has blossomed in a large variety 
of physical systems, from atoms [2] to semiconductors [3] 
passing by superconducting circuits [T and nanomechan- 
ical oscillators [5j. A fascinating aspect of this funda- 
mental problem is how it bridges the gap between quan- 
tum and classical coherence. In the former case, one has 
quantum superpositions of light and matter, entangling 
photons with the ground and excited states of the emit- 
ter. In the latter case, one has a classical photon field, a 
continuous function of a continuous variable, fully spec- 
ified in all its attributes. The passage from one to the 
other can be tracked in the one-atom lasing transition. 
The one-atom laser is a concept first proposed and the- 
oretically studied by Mu and Savage [6] , with the aim of 
achieving lower thresholds for lasing. They encouraged 
experimentalists to bring the number of atoms in a 
conventional laser to unity. In a very high quality fac- 
tor cavity, the emitter reaches the strong coupling regime 
at the single excitation level. They showed that in this 
regime, a single incoherently excited emitter (a two level 
system in the simplest case) can constitute the whole gain 
medium and populate singlehandedly the cavity with a 
very large number of photons. If the spontaneous emis- 
sion rate of the atom into other modes than the cavity is 
small, the growth in the population of photons exhibits 
no threshold as a function of the pumping rate and 
develops a classical (Poissonian) field statistics, thanks 
to the efficient periodic exchange of excitations with the 
atom. 

On the theoretical side, the single-atom laser has been 
extensively studied for its qualities as a laser [SI |H1- 
[22], mostly solving its steady state numerically but 
also through analytical techniques, such as the contin- 
ued fraction expansion |22j or phase-space representa- 



tion [121 HSl [22], and under different approximations 
such as the few-photon [21 dH 1201 112] or semiclassi- 
cal [S] [ni dZl mi US dynamics. On the other hand, 
the spectral properties, that is, atom and cavity photo- 
luminescence emission spectra, have been studied only 
numerically [ilHllSl HH Ullll dl] , although some ap- 
proximations for the linewidth of the cavity spectrum, 
that converges to a single peak and exhibits the stan- 
dard line-narrowing, provided useful analytical expres- 
sions [m |T7]. Given that the cavity field is coherent 
in the lasing regime, the back-action of the field on the 
emitter leads to the formation of a Mollow triplet in its 
spontaneous emission spectrum [9] [Ml [16] , with similar 
properties to that theoretically predicted by Mollow for 
an atom under resonant laser excitation [23] . This struc- 
ture is of a great fundamental interest as it represents a 
pinnacle of nonlinear quantum optics. 

On the experimental side, the strong coupling regime 
is now firmly established at the single and few photon 
level with atoms [24' as well as with artificial atoms, su- 
perconducting qubits ^25. ,26 or semiconductor quantum 
dots [27H29] . and the one atom laser as described above 
has been realised in all these systems [501132] . The Mol- 
low triplet under incoherent pumping has not yet been 
reported in any experiment, one reason being that they 
typically focus on the cavity field lasing properties. Only 
the original configuration proposed by Mollow, under res- 
onant coherent excitation, has been reported, also in all 
the systems above, namely with a single atom ^33,, a 
molecule [SJj ^ superconducting qubit [35] and a quan- 
tum dot [M3S- 

In this text, we consider in detail the quantum to clas- 
sical transition that leads to lasing in strong-coupling, 
starting from the fully quantized description of Jaynes 
and Cummings [40j . We show that quantization breaks 
down as coherence is formed and a classical description 
becomes more appropriate. We provide several limiting 
cases that describe the system in its various regimes of 
operation. Varying degrees of agreement are afforded de- 
pending on the complexity of the approximation. We 



2 



provide compact analytical approximations in the sim- 
plest cases and a straightforward numerical procedure 
that leads to an excellent quantitative agreement. We 
include pure dephasing, important in solid state sys- 
tems [1TJI12], and arbitrary detuning between the modes. 
We focus more particularly on the lasing regime at res- 
onance, where we show that the emergence of a MoUow 
triplet manifests classical nonlinearities of strong light- 
matter coupling in cavity-QED. As a whole, we show 
that the lasing transition is a rich and complex one and 
we hope to give a rather comprehensive view of its various 
limits. 

The remaining of this text is organized as follows. As 
the most striking and characteristic manifestation of las- 
ing in strong coupling is the MoUow triplet formed un- 
der incoherent pumping, we first revisit its coherent ex- 
citation counterpart, this is done in Sec. |Tlj extending it 
to include pure dephasing and detuning. Starting from 
Sec. |III[ we turn to the case of incoherent excitation ex- 
clusively and obtain analytically the steady state, mode 
populations and photon counting statistics ( III A ) , the 



system full density matrix ( III B I and the two-time cor- 



relators needed to compute the power (or luminescence) 
spectra (IIIC). We derive the expressions for the lasing 



properties by applying the semiclassical approximation 
that we compare to other approximations that describe 
the transition into and out-of lasing. In Sec. |IV( we put 
all these elements together and derive the spectra of emis- 
sion for both the cavity and the emitter. We analyse the 
resonances of the system ( IV B ) as well as the elastic scat- 



tering component (IVC). We apply again the semiclas- 



sical approximation to simplify the emitter spectra into 
a compact closed-form expression for the MoUow triplet 
under incoherent pumping ( IV D ) . This expression is used 



to explore the parameters where the MoUow can be ob- 
served experimentally. In Sec.|Vj we summarize our main 
findings. 



can be removed by going into a frame rotating with the 
laser (A = wl — Wo-): 



(2) 



The spontaneous decay and the pure dephasing suffered 
by the emitter can be described by two Lindblad terms 
in the master equation: 



dtp = i[p,Hi,] + 



Ha 



70 



(3) 



where Ca{p) — {2ap(7^ — a^ap — pa^a). The steady state 
of this simple system {dtp — 0) can be solved analytically 



(see Appendix 
coherence ESI 



A|) in terms of the emitter population and 
431: 
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The effective coupling to the laser (intensity that effec- 
tively excites the emitter) is: 
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At resonance, — flh and (u^) is pure imaginary. The 
effective laser intensity fijf is reduced with detuning by 
an amount that depends on the overlap in frequency be- 
tween the laser and the emitter lineshapes. As the laser 
has no linewidth, the total emitter width 7^ -I- 7^ deter- 
mines the overlap. Therefore pure dephasing compen- 
sates for detuning by increasing this overlap. 

The normalized spectra of emission reads in the steady 
state (that we set as i = 0): 
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II. MOLLOW TRIPLET UNDER COHERENT 
EXCITATION 

The analysis of light-scattering by a two-level system 
(representing an atomic transition) was first given by 
MoUow [23j , who reported the antibunching of the scat- 
tered light, as well as the spectral structure now known 
as the MoUow triplet [33]. It results in the case where a 
strong-beam of light, with frequency wl, impinges on the 
emitter which has a natural frequency uj^ (^ = !)• The 
Hamiltonian reads: 



(1) 



where a is the pseudo spin operator for the two-level 
system and flh is its coupling strength with the optical 
laser field. Note that the later, described by a complex (c- 



number) wave E{t) w f^L(e 



), is thus entirely 



classical. The explicit (and fast) time dependence in Hj^ 



Any two-time correlators can always be decomposed 
as a sum of complex damped exponentials |44j : 



(at(O)a(T)) = 



p 



(7) 



where all the parameters, weights Lp, Kp, frequencies 
ujp and effective decay rates jp, are real. They can be 
obtained by means of the quantum regression formula 
(see Appendix [A|). Eq. ^ leads to the spectrum, 
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In the case of Eqs. (ll||3|), the spectrum has four compo- 
nents, that we label p = coh, 0, +, — . The elastic scat- 
tering component (p = coh) is a delta peak 6{uj) at the 
laser frequency (ojcoh ~ 7coh ~ Kcoh = 0) with weight: 
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FIG. 1: (Color online) (a) Origin of the peaks in the MoUow 
triplet: the three different frequencies (p = 0, ±) are found in 
the four possible transitions between two Jaynes-Cummings 
rungs at high intensities. Below, different MoUow triplets 
when varying detuning (b, d) or dephasing (c), with JIl = 
1.57cr. The symmetry of the triplet is broken only under the 
combined action of detuning and dephasing. 



Strong coupling, where the character of the dynam- 
ics of the two-time correlator is oscillating rather than 
damped, is defined by the appearance of this splitting 
(3?(i?L) ^ 0), that is: 



2r!L > |7^- 701/4. 



(12) 



We see from Eqs. (10b I and (lOd) that, beyond the 



expected broadening of the lines, dephasing also shifts 
the two satellite peaks. In general, this shift brings the 
side peaks closer to each other, inducing the transition 
into weak coupling when 7,/, > 70- + SJIl . However, sur- 
prisingly, the maximum splitting, for a fixed 7a- and JIl, 
corresponds to a nonzero dephasing, 7^ = 7^. In fact, 
the splitting remains different from zero, in the presence 
of dephasing, as long as 7^ — SJIl < 7^ < ^a + '&VL'L. If the 
driving field is too weak to bring by itself the system to 
strong coupling (JIl < 70- /8), it can be aided by increas- 
ing dephasing. However, higher dephasing also blurs the 
spectral features. For the regimes of excitation where de- 
phasing induces strong-coupling (0 < 7^ — 8f2L < 7^ < 
7ct), the observed lineshape always remains single peaked. 

The final expression for the MoUow triplet spectrum 
at resonance in presence of dephasing reads: 
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The inelastic scattering part is a triplet with a central 
peak {p ~ Q) and two sidebands p = ± that carry the in- 
formation about the light-matter interaction. The physi- 
cal origin of these two peaks is in the transitions between 
the dressed states of the Jaynes-Cummings Hamiltonian 
at high number of excitations IS' , as shown in Fig. [ija) . 
The two transitions between different types of dressed 
states become degenerate and form the central peak (0), 
while transitions between the same type of dressed states 
give rise to the side bands. The full expression of the 
spectrum out of resonance is too lengthy to be given here. 
The resonant case formula is shorter. It reads: 



in 
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t^o = , 70 = 7ff + 70 : 
L± + iK± = 



(10a) 

(10b) 
(10c) 
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where we have defined the (half) Mollow splitting: 
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The S scattering peak and the central peak in the first 
line are neatly set apart from the two side bands in the 
rest of the expression. This decomposition is shown in 
Fig. Ha). 

The effect on the Mollow triplet of detuning the laser 
from the emitter is shown in Fig.[ljb). It spreads the side 
bands apart, with asymptotes tOa- and Wo- -I- 2A, while the 
central peak, pinned at the driving laser frequency wl, 
gets suppressed. The scattering peak, not shown, ulti- 
mately dominates the spectrum over its incoherent part 
which fades away. In any case, the lineshape remains al- 
ways symmetric with respect to the laser frequency (cen- 
tral peak), a characteristic proper to coherent excitation 
as we shall see later. 

Out of resonance, pure dephasing has a strong qual- 
itative effect: it breaks the above symmetry, bringing 
the spectrum towards the uncoupled case with only one 
Lorentzian peak at Wo- with FWHM 7^ +70- Even a 
small dephasing enhances considerably the emitter peak 
relatively to the others, as shown in Fig.[l|c). This asym- 
metry becomes larger the weaker the effective laser drive, 
that is, for lower f^L and larger A and dephasing. One 



4 




FIG. 2: (Color online) Degree of asymmetry V of the MoUow 
triplet under coherent excitation, as a function of detuning A 
and pure dephasing 7^ for JIl = l-5ya- It is (dark blue) 
when the side peaks have the same intensity and 1 (bright 
yellow) when one of them disappears completely. 



can quantify the visibility of this asymmetry as the dif- 
ference between the intensities of the two side peaks: 



V = 



\L+-L_\ 
\L,\ + \L^ 



(14) 



This is plotted in Fig. [2j where lighter (yellow) colors 
refer to smaller degrees of symmetry (minimum when 
only one peak of the two side bands survives). As a 
practical application, one can measure the magnitude of 
pure dephasing as a function of detuning from the degree 
of asymmetry. 



III. ONE-ATOM LASER 

In the cavity QED version of this physics, the system 
is described by the Jaynes-Cummings Hamiltonian [45]: 



(15) 



where also the light field is quantized, through the anni- 
hilation operator a. The detuning is now A = coa — uJa- 
and we consider = as the reference energy. The 
Liouvillian, dtp = Cp, to describe this system in a dissi- 
pative context with decay (7c), incoherent pumping (P^) 
and pure dephasing (7^) has the form [T5]: 

£0^i[0,H]+ ^{2cOc^ -c^cO-Oc^c) (16a) 

c—a,a 

+ J2 ^{'2c^Oc-cc'0-Occ<) (16b) 



(16c) 



where p is the density matrix for the combined emit- 
ter/cavity system. The effective broadenings of the un- 
coupled modes are defined by Fj, = 7^ — Pa and F^ = 



A. One-time correlators: populations and statistics 

The light field that was previously a classical laser field 
was fully characterised by its intensity (|riLp) and its 
frequency (wl). In the fully quantized description, corre- 
lations between the fields should be taken into account, 
namely, in the steady state: 



a"- W) , (17a) 



(17b) 



with f^a and N„ \n\ real and iNaa \n\ complex in general 
but real at resonance, all others being zero. The main 
observables that characterize the system are: 

= 7Va[l], n,^N,\\\ and g'^) ^ iV„[2]/n2 . 

(18) 

In the following, we provide exact implicit expressions 
for the correlators (17), that allow an efficient numerical 



solution, and derive approximate analytical expressions 
for different regimes of excitation. 

In the case without any direct cavity pumping, = 0, 
the field correlators admit a simple expression in terms of 
Na\n\ (the general equations are given in Appendix [B|) : 



+ 7a(«- 1) 



2.9 



Tcr + 70 +7a(2n - 1) ' 

This allows to obtain a single equation for Na \n\ 
1 , n-ia 2P, 
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To- + (n - l)7a 



F^ + (n - 1)7q 



(19a) 
(19b) 

(19c) 
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Na[n + 1]. 



(20) 



where we have introduced, respectively, the effective co- 
operativity, effective coupling and total decoherence (in 
the presence of detuning and pure dephasing): 



v/1 + (2A/FtN)2' 
FtM =Ta + 70 + (2?^ - l)7a • 



(21a) 

(21b) 
(21c) 



Let us note that ^off [n] = g for all n at resonance or when 
decoherence is large as compared to the detuning. As in 
the case of laser excitation, the effect of detuning is to 
effectively diminish the coherent coupling, which mag- 
nitude is linked to the spectral overlap between modes 
(represented by Fx). Here too, the decoupling caused by 
detuning can be compensated by increasing decoherence 
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FIG. 3: (Color online) Exact results computed numerically for (a) ^ana/g, (b) and (c) g'^', as a function of pumping 
for various systems, with 7a/g G {0.1,0.22,0.46, 1,2.15,4.64, 10} from lighter to darker shades (top to bottom in (a)). Other 
parameters are 7ct = 0.00334g, from Ref. |47], and Pa = 7</> = A = 0. The different regimes of operation are designated in (a). 



(decay or pure dephasing) since in this case the spec- 
tral overlap between the cavity and the emitter increases, 
bringing them effectively back to resonance. Detuning 
and pure dephasing only appear in the cooperativity pa- 
rameter Coff [n] , as noted also by Auffeves et al. [20] . This 
situation is similar to the case of two coupled harmonic 
modes [iS] . 

We reduced the whole steady-state problem of the 
Jaynes-Cummings with emitter pumping and decay to 
a single equation, Eq. ( 20 1 , which is however a nonlin- 



ear recurrence equation with non-constant coefficients, 
for which there is no general method leading to an exact 
solution. The problem put in this form is nevertheless 
quite tractable numerically and we shall in the follow- 
ing present various limiting cases which will spell out the 
physics of this problem. 



Since Na[0] = 1 by definition, Eq. (20) can be easily 
iterated numerically to provide Na [n] for all n as a func- 
tion only of the mean number of photons in the cavity, 
Ua- The small n equations are the most important since 
they capture the dominant few-photons correlations. The 
exact expressions for and in terms of Ua arc, in any 
case, simple enough: 



i(2) - 



Pg - Ian a 
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(22a) 
(22b) 



2P„ 



+ 70 + 7a 7a + ■ 



7a 



They are given in terms of a key parameter of the system, 
the Purcell rate of transfer of population from emitter to 
the cavity mode: 



4(gcff[l] 

7a 



(23) 



This parameter is large for good cavities, when cavity 
QED is realized at its fullest: k„ ^ g. 

One can obtain Ua self-consistently by truncating 
Na[n] at a sufficiently high number of photons, rtmax, and 
solving numerically the resulting finite set of equations. 
This gives the results plotted in Fig. [s] for (a) 7a'T-a/3, 



(b) n„ and (c) g^'^\ In the best systems (7a,7cr <C g), 
we distinguish five regions in these plots, which we shall 
investigate in more details in the remaining of this text: 

1. Linear quantum regime (or simply ^Hinear''''), where 
Per <C 7ct, keeping the emitter essentially in its 
ground state, very rarely excited. 

2. Nonlinear quantum regime (or simply quantum^^ ) , 
where ^ is enough to probe higher (n > 
1) rungs of the Jaynes-Cummings ladder, without 
climbing it too highly so that few photons effects 
remain the dominant ones. 

3. basing or nonlinear classical regime, when 3> 7^ 
and the emitter population w 0.5, the cavity can 
accumulate a great number of photons and the field 
becomes Poissonian. 

4. Self- quenching regime, when P„ ^ K(j/2 starts to 
drive the emitter to saturation, Ua > 1/2, reducing 
the number of photons. 

5. Thermal regime or linear classical regime, when 
P^ > K„, the emitter is always in its excited state, 
rio- — > 1, and the dephasing induced by the pump 
disrupts the coherent coupling, so that the number 
of photons is very low again and the field becomes 
thermal. 

Similar classifications have been proposed, for instance 
by Poddubny et al. [l7]. In the following, we address 
several types of approximations, that perform to varying 
degrees of accuracy depending on the level of complexity 
involved and the regime under consideration. Some of our 
approximations recover known results [HI HU [T7| l^ll HI] , 
however, this comparative analysis will give us, beyond 
good approximated formulas, valuable insights into the 
underlying physics. It will also allow us to determine 
the pumping ranges that determine each regime. As the 
expression for n„ and g*-^' follow straightforwardly from 
that of na , we will not provide their explicit form in most 
of the cases analysed below. 



6 



(a) 



/ X, v'-,7-^- - -\ 1. \, \ \ 

.-■'1 ■f' .y '^'^f-' i-C";- •• 



boson (linear model) 
truncated JC 



(b) 




(c) 



thermal field approximation 




0.001 



1000 



Pa/g 



FIG. 4: (Color online) Comparison between the numerical 
results for Ua (from Fig. [sja)), in thick yellowish lines, with 
different approximated solutions. In each panel, 7^ decreases 
exponentially from bottom to top curves, ranking from lOg 
(darkest curve) to O.lg (lightest curve). The limiting case 
7a = is also shown on the left upper corner, where the di- 
vergence is a feature of the Jaynes-Cummings model. In (a), 
we superimpose the linear model (dotted red) and Jaynes- 
Cummings truncated at one excitation (dashed blue), given 
by Eq. ( |24||(a ). In (b), the approximated semiclassical solu- 
tion. Eg. ( |32[ |, which provides an accurate description in the 
lasing regime given by Eq. (61 1. In (c), the thermal approx- 
imation, given by Eq. ( 39 1 that converges to the models in 

(d), _ 

In this 



(a) at low and high pumps. In (d), the cothermal approxi- 
mation, given by the numerical solution of Eqs. 



(441 



case, we extend solutions to values of "/a/g till O.Olg, out of 
reach numerically. A vertical guideline marks the value of 
7ct = 0.00334(?. Other parameters are = A = 7^ = 0. 



1. Linear model approximations 

In the linear regime, where the emitter is excited with 
very low probability, it can be well approximated by 
another harmonic oscillator [3S]. This allows to find a 
closed-form analytical solution: 



l^aiTa + 7a) + ^aiTa + 7a + 1^) 
K<T + 7a + 7(T + 70 



Pa, 



(24a) 
(24b) 



Equivalent expressions for and are obtained in 
the first order truncation of a continuous fraction ex- 
pansion of these quantities, as recently shown by Gart- 
ner [22j. They are also formally identical to those ob- 
tained by truncating the Jaynes-Cummings model at the 
first rung of excitation [201 113 HH] [Si] . The only differ- 
ence is in the effective broadening T„, that appears with 
a — sign in the case of coupled bosons, 7^ — F^. (48j . 

and a -I- sign in the truncated Jaynes-Cummings model, 
Fcr —> 7o- + Per- At Pa <C 7ct, the sign becomes irrelevant 
with Fo- w 7o- and the population grows linearly with 
pumping, Ua ~ CiPa, with the slope: 

Ci = ^ • (25) 

Kcr(7<T + 7a) + laKla + 7a + 70) 

This agrees with the numerical results for P^ ~ 7cr as 
shown in Fig.ga). In the case where 70- = 0, we simply 
have Ua ~ Pa Ma- Interestingly, the two models also 
provide the right formula in the high pumping regime 
where the number of photons is low again, 



lim na = — 

P„^OD P„ 



(26) 



but the emitter is completely saturated. 

The bosonic populations diverge at two values of 
pumping (where the denominators vanish): 



P±=la 



■ 7a + 70 



i± Wi- 



4KCT7a 



Kct + 7a + 70 



(27) 

For good systems with small cavity decay rates, we have 
P_ ~ 7o- and P+ w + 7^ + 7^. At P_, the two 
populations diverge but remain positive, a manifestation 
in this model that the system enters the lasing oscil- 
lations where both populations are "inverted" . For in- 
termediate pumpings, where the system is in the lasing 
regime P- < P^ < P+, both populations are negative, 
meaning that the physics in this regime is out of reach 
of this model. For large enough pumping, P^ > P+, 
the system exits the lasing regime as in the bosonic 
model converges again to the exact numerical result fol- 
lowing Eq. (26 1. However, ria remains negative for all 
Pct > P-, meaning that its population remains there- 
after inverted. The bosonic model thus provides an ac- 
curate description of the transition in and out of lasing 
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though the appearance of divergences and negative pop- 
ulations, as we will confirm later when linking it to the 
full Jaynes-Cummings system. The truncated Jaynes- 
Cummings formulas remain always positive and instead 
of a divergence, the two-level system becomes inverted, 
with — 7- 1 from below. Given that the number of 
photons is always below one, this model is not suited to 
describe the transitions in and out the lasing regime. 

The bosonic and fermionic models provide opposite 
statistics for the cavity field. Two coupled harmonic os- 
cillators under incoherent excitation are always thermal 
with a photon distribution T{n) = ii?'^/{na + 1)"+^ and 
(^(2) _ 2. On the other hand, truncating the Jaynes- 
Cummings model at one excitation means that one ex- 
cludes the possibility to have two photons at a time in 
the cavity, which results in perfect photon antibunching. 
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(2) 



0. Both models provide the exact Ua solution of 



the Jaynes-Cummings model in two and opposite limit- 
ing cases, namely, 7a for the bosonic model, where 



(28) 



la -Pa' 



with g(^) = 2, and 7a — > oo for the truncated Jaynes- 
Cummings model, where 



(29) 



la + Pa la' 



with f/(^^ = 0. This is seen in Fig. |4|a), where 7a/<? = 10 
(lowest curve) is already a good approximation for cxd, 
while 7a/5 = (uppermost curve) is exact. In the first 
limiting case, recovered by the bosonic model, the system 
accumulates photons in the most effective way possible, 
which leads to a divergence in the number of photons, 
also in the exact Jaynes-Cummings model. The Rabi 
delivery of photons is so efficient in the one-atom laser 
that, unless there is a leakage of photons from another 
channel, the accumulation of photons is unbounded. It 
is unlimited by the strong-coupling feed. 

In the opposite limit of weak and inefficient coupling, 
the emitter undergoes population inversion unaffected 
by the cavity while the cavity gets an effective pump- 
ing of photons from n„ through a very weak Purcell rate 
K(j [20 , 22 . For 7a > 4^ the truncated Jaynes-Cummings 
model provides a quantitatively good agreement for the 
entire range of pumping, as shown in Fig. [4][a), if one 
excludes the behaviour of g^'^\ This limit is also ac- 
counted for exactly by a series expansion in P^, given in 
Appendix [Cj since the coupling is perturbative. In this 
case, the system goes directly from the quantum linear 
to the classical linear regime. 

One can obtain an exact expression for g^^^ in the linear 
regime solving Eqs. (20) truncated, not at one, but at two 
photons, that is, for n = 1, 2 (assuming iVa[3] = 0). Note 
that at Per = 0, Na[n\ — for all n, but we are interested 
in the limit: 



(2) 
dpi- 



lim 



Na[2] 



(30) 



that remains different from zero, since Na[n] cx P". This 
method leads, to first order in P^, to the populations of 



Eq. (24) and to: 



(2) 



Kaila + la) + 7o-(7a + la + l4>) 



l^ai^la + la) + (7a + la){ila + la + I4,) 



(31) 

This result is exact and valid for any set of parameters, 
as one can check by simply solving the equations to the 
next order of truncation (n = 1,2,3 and iVa[4] =0). In 
general, g'-^-' £ [0, 2] in the linear regime of the Jaynes- 
Cummings dynamics, although in Fig. ^c) we only show 
< s''^-' < I since we have chosen 70- ~ 0. 



2. Semiclassical approximation 

Given that, as shown in Fig. |3][c), the cavity field be- 
comes Poissonian in the lasing regime, we can find good 
approximated solutions for the steady state under the as- 
sumption T[n] = e~'^'^n'^/n\, which leads to A^aM = n^- 
To establish the lasing in strong-coupling regime, one also 
needs a good cavity, so we shall assume 7a ^ 5, and high 
enough pumping, Po- » 7^, 7a- Plugging Na[n\ = rt" into 
Eq. (20) and solving the resulting equ ation for n = 1, or 
equivalently, imposing g^^^ = 1 in Eq. (22b), gives 



27a 



1 



27. 



■70 



(32) 



This is a very good approximation for the region where 
the cavity field behaves classically [HI |T2] , equivalent to 
solving the n = equation in Eqs. ^ ^ [52]. The 
probability of finding the emitter in its excited state 
reads: 



1 



1 



-70 



(33) 



The approximated ria is plotted with blue dashed lines 
m Fig. igb) for comparison with the numerical results 
and exhibits a remarkable agreement for a large pump- 
ing range of practical interest (the axis is in log scale). 
Similar agreement is found for |21| . 

The two expressions for the populations have a 
straightforward interpretation. Two parameters deter- 
mine the populations (neglecting for the sake of sim- 
plicity the small correction brought by the emitter de- 
cay, 2^cr/^a ^ 1): the cavity feeding''^ and the emit- 
ter feeding^^ efficiencies, defined as Fa — T^/{2^a) and 
Fa = (Ta +l<ji)li^a, respectively. They follow as: 



Fa{l~F,) and «(l + P,)/2, 



(34) 



The cavity population increases linearly with pumping 
[ua ~ Fa) while the emitter is half occupied (rig. « 1/2). 
This is the range of pumping with the most effective ac- 
cumulation of photons in the cavity, as the incoherent 
processes are small enough not to disrupt the coherent 
coupling dynamics. All the excitations injected into the 
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emitter are transferred into the cavity. We already have 
seen that there is a hnear relationship between Ua and 
Pct hi the — aptly denominated — linear regime, as given 
by Eq. ([25]). A similar linear relationship Ua ~ 6*2^0- 
also holds in the lasing regime, when Pa- but be- 

yond the quantum regime, > ^aila, this time with a 
slope C2 defined as: 



Co = 



1 

27^' 



(35) 



The transition between the two types of linear behaviours 
ria ~ CiPa, i = 1 — >■ 2, and the question of the threshold 
in this process is an interesting subject |14j that would 
bring us too far astray and that we postpone to another 
work [SO] • We will only comment in the present text that 
this intermediate region is the less liable to the types of 
approximation that we derive here, since it lies at the 
frontier between the very few and the very large num- 
ber of excitations, and no good approximation can re- 
produce it to a high degree of accuracy other than by 
keeping track of all correlations between the particles, 
which, being a N body type of problem, implies a nu- 
merical procedure. We call this intermediate region the 
"quantum nonlinear" or simply the "quantum" regime. 

When the pumping is sizable as compared to k^-, the 
emitter occupation starts to show signs of saturation, in- 
creasing linearly with pumping, and quenching the linear 
increase of the cavity population. represents therefore 
the degree to which the pumping succeeds in populat- 
ing the emitter itself, against the coherent exchange of 
population that feeds the cavity with efficiency Fa- The 
maximum population of the cavity 



max(na) « ^(1 " ) 

"7o l^a 



is reached at the intermediate rate 

27. 



\ max.{na ) 



1- 



70 ^ 



(36) 



(37) 



The present approximated expressions are valid until no- 
approaches 1, then, the self-quenching dominates the dy- 
namics and ria — >■ 0, at around: 



Pr, 



Sjcr - 270 , 



(38) 



that is, when the pump reaches the effective transfer rate 
of excitation towards the cavity mode. At this point the 
statistics changes to thermal due to pump-induced deco- 
herence. Note that Pmax ~ P+, found in the previous 
subsection when analysing the bosonic model. 



3. Thermal approximation 

Assuming a thermal state for the cavity field, with 
statistics T{n) = nJlliua + 1)""*"^ and thus with Na\n\ = 
n!n", satisfies in good approximation Eqs. (201 when the 



system undergoes self-quenching and gets driven into the 
thermal regime (P^ > k^). The Ua obtained from the 
equation (n — 1), or equivalently, imposing g^^-* = 2 in 
Eq. (l22bl), reads: 



-{(r. + 7a)x 



87a 



2P„ 



16Pa7a/ra ^ n ^ + 7a + 70 ^ 

Tff + 7a Ha Va + 7a 

-r.(^:^ + ^-i)-7a(^^^ + 2) 



7a 1 2 



r 



(^ + 1)}. (39) 



This solution converges to the linear result of Eq. (24) 
to first order in P^ when Pa < ^a- The opposite limit, 
well into the thermal region, is that already found with 
Eq. (26) to first order in 1/Po-, which provides the de- 
creasing tail at large pumpings, as shown in Fig. |4j The 
lines converge to the same universal curve when plotting 
la^al 9 in Fig. [3]ja). In the intermediate region, it pro- 
vides a good approximation when the system is not good 
enough to lase, at large dissipation rates, similarly to the 
truncated Jaynes-Cummings formula. Of course, having 
— 2 throughout, does not, in general, represent well 
the statistics. 

Fig. [ijjc) may seem to shown that Eq. ( 39 ) captures 
the qualitative behaviour of Ua even in the lasing regime. 
It does not, however, provide the change into the linear 
lasing slope from ^ C\Pa to C2Pa, discussed in the 
previous subsection, which is a serious conceptual short- 
coming. 

In two particular cases, the thermal solution becomes 
exact for Eq. (20), namely, 7c = Pc = with c — a (pre- 



viously discussed) on the one hand and its counterpart 
c = cr on the other hand. In both cases Naa = and 
Na[n] = fiaNaln - 1] with Na[n] = n\nl 



''n", where: 



Pr 



Ic-Pc 



and 



Ic+Pc 



(40) 



The steady state in this case is two uncorrelated ther- 
mal fields at the same temperature. This is independent 
of the coupling strength g, which only determines the 
speed at which this steady state is achieved. Although 
the discontinuity at 5 = might seem counterintuitive, it 
is physically clear that the thermal equilibrium does not 
otherwise depend on details of the microscopic couplings. 
The case Po- = 7^ = 0, that corresponds to thermal ex- 
citation of the cavity mode only, presents interesting as- 
pects which we also postpone to a future work. 

The case Pq = = 7a = 0, appearing in Fig. [4]jc) as 
the uppermost exact curve, diverges, as we noted pre- 
viously, at Pa > 7cr, meaning that the system does not 
have a steady state but rather its intensity grows with- 
out bounds. This is to be interpreted as the instability 
accompanying the transition into lasing [12] . 
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1000 



FIG. 5: (Color online) Comparison between the numerical 
g'^' (in thick lines) and the cothermal approximation (dashed 



lines), from Eq. (43 1 with Ua of Fig. |4] It interpolates neatly 
between the lasing and the thermal regions (1 — >■ 2) but also 
provides a very good agreement in the linear regime. 



4. Cothermal approximation 

A very good quantitative approximation for the whole 
pumping range is obtained assuming that the field is in a 
so-called cothermal state [51 , which is the field with, in 
average, rf°^ coherent photons and n*'* thermal photons 
(with a total cavity population na = n^a + "a°^)- this 
case, the photonic statistics reads: 



T[n] 



-n, 



coh 



(l + <'')»+i "^<^(l+n*^)^' 

(41) 

and moments of the distribution are: 

^coh 



-ith 



(42) 



with L„ the n-th Laguerre polynomial. The second order 
coherence function is 
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(2) = 2 - 



(43) 



One can write the set of Eqs. ( |20[ ) for n =1, 2 using the 
parametrization of Eq. (42): 







1 



Ccff[l] 



Pa 



27a 



7a' 



2P. 



1 



-7a 



= 



1 



27a 



2P. 



■Ccff [2] 
2P. 



7a 



f 27a 

27a 
. +27a 



1 



(44a) 
[2r.a - «°')'] 



+ 27a 

[6n3_9na«°'>)2+4«°'>)3]. 

(44b) 



Equations (44) can be easily solved numerically (two 
coupled nonlinear equations) , extracting the physical so- 
lution Ua, njj"^ > 0. They yield an excellent agree- 
ment with the complete numerical solution as shown in 



Fig. ^d) for the cavity population and Fig. 

[5]for5(2). We 

took advantage of this to extend the results to cases out 
of reach of our numerical procedure, namely, ja/g down 
to 0.01, shown on Figs. |4] and [5] as the dashed lines only, 
without the corresponding numerical calculation. The 
transitions from the linear regime to lasing and from las- 
ing to the thermal regime are well accounted for within 
this approximation. Although they are not exact, they 
provide a fairly good quantitative description of the en- 
tire pumping range for any set of parameters, particularly 
of the quantum regime which eludes most approximation 
schemes. 

The transition into a thermal field can be induced by 
increasing the incoherent pumping but also other sources 
of decoherence such as dissipation, pure dephasing or de- 
tuning. In Fig. [g] we plot the extracted ria and g'-^^ 
respectively, under the cothermal approximation, as a 
function of pumping, when increasing detuning and pure 
dephasing. The lasing to thermal transition is even 
more apparent in the Mandel-Q parameter, defined as 
Q ~ "alff^^-* — 1), featuring it with a broad peak, as shown 
in the third row of Fig. [6j Both detuning and dephasing 
are clearly detrimental for the intensity and coherence in 
the self-quenching regime but they may compensate each 
other in the transition from the linear regime into lasing. 
Fig.[6]shows the appearance of a threshold in the presence 
of detuning (at small pumping the thermal fraction in- 
creases) , due to the reduction of ^eff [Z] ■ Adding dephas- 
ing, 70, may compensate this reduction and smoothen 
the threshold up to a certain point (as it also increases 
to the total decoherence rate Ft). This observation was 
made by Auffeves et al. |20j . In general, the maximum 
intensity, max(na), and coherence achieved, always de- 
crease with both detuning and dephasing. However, de- 
tuning makes it occur at higher pumpings while dephas- 
ing does at lower ones. Pure dephasing is thus indeed 
making both transitions, from quantum to lasing and 
from lasing to self-quenching, occur at smaller pumping. 



B. Density matrix 

We have obtained in the previous section good approx- 
imations for the statistics of the cavity field, T[n], as well 
as the main quantities of interest such as Ua- To compute 
the photoluminescence spectrum, what we will do in next 
section, we also need the full density matrix in the steady 
state. In the very strong-coupling regime, we can relate 
it analytically to T[n], as we show in this section. 

The full statistics is most conveniently obtained from 
the master equation with elements Pm,i;n,j for m, n 
photons and z, j excitation of the emitter (m, n G N, 
i,i E {0, 1}). Rather than to consider the equations of 
motion for the matrix elements directly, it is clearer and 
more efficient to consider only elements that are nonzero 
in the steady state. These are: 



Pri,0:n,0 i Pl["-] — Pn,l;n,l , — Pnfi;n~ 



-1,1 1 

(45) 
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A/^ = 4 A/^ = 8 A/^ = 12 




Pa/g Pa/g Pa/g Pa/g 

FIG. 6: (Color online) (a) Ha, (b) (;^^' and (c) Q-Mandel parameter, as a function of the pumping rate, extracted with the 
cothermal approximation. Detuning is fixed in each column to the value while pure dephasing increases from (thick curves) 
to 12g (dashed curves) in steps of 2g. Detuning brings a threshold for lasing and while dephasing can help to reduce it, it also 
makes the system less robust to quenching. Parameters are 7a — O.lp, Pa = = 7o- = 0. 



and correspond to, respectively, the probability to have n 
photons with (pi) or without (po) excitation of the emit- 
ter, and the coherence element between the states \n, 0) 
and |n— 1,1), linked by the Hamiltonian Eq. (15 1. 
Both Po and pi are real. It is convenient to separate q 
into its real and imaginary parts, q[n\ = (ir[n\ + i<li[n\, 
as they play different roles in the dynamics. The equa- 
tions for these quantities are derived from the Liouvillian 
equation (16). They read (their full form is given in Ap- 
pendix Ib|) : 



dtpoin + 1] 

dtPi N 
dtqi[n + l] 



=2?phot{po[n- + 1]} 

- PcrPo[n + 1] +7crPl[^ + 1] 

- 2gVnTlqi[n + 1] , 
=2?phot{Pi W} 



- 7o-PiN + PaPo[n] - 

=^?phot{qi["- + 1]} 

T,+7^ , 2A2 



(46b) 

Tqi[n+ 1] , 



■70 



+ gVn + l{po[n + 1] 



qi[n + l] 
PlW), 



(46c) 



where we have separated the cavity dynamics into a su- 
peroperator Pphotj which exact expression is given in the 
Appendix. If this photonic dynamics is much slower than 
the rest, that is ja,Pa ^ ^<7,g, one can solve the emit- 
ter dynamics separately, ignoring Pphot El] ■ This as- 
sumes that the photon distribution T[n] = poM + Pi[n\ 



does not change during the excitation and interaction 
with the emitter that happens at a timescale, of order 
1/Pa or 1/g, much faster than the photon one, of order 
l/7a. This approximation becomes exact for the per- 
fect cavity, when 7a = (and Pa = 0, which we already 
assumed), which is why the system admits an exact so- 
lution (cf. Eq. (40 1). Neglecting the photon dynamics 



allows one to link all the density matrix elements with 
the photon statistics T[7i]: 

Po[n + l]w (47a) 



(46a) ^^(„ + i)(^^T[n] + ^T[n + l]) +7,T[7i + l] 



2K„(n + l)+r, 
PiN~ (47b) 
t^ain + 1) (^T[n] + ^T[n + 1]) + P,T[n] 



qi[n + 1] 



2Ka(n + l)+r^ 

KaVn+ 1 



25 



(po[n + l] -pi[n]) (47c) 



where 



KaVTr+T P,T[n] - jaT[n + 1] 
2g 2Ka{n + l) + r^ 

4(.9eff[l])2 



p ^ (48) 
Tct + 70 

is the Purcell rate of transfer of population from the cav- 



ity mode to the emitter (cf. Eq. (23l) and gos follows 
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the definition in Eq. (21 1 for negUgible 7^. Note that 
Po[n] + Pi[n] is not strictly equal to T[n], due to our ap- 
proximations, but the numerical discrepancy is small in 
the regime of interest where the number of photons is 
high and n ~ n + 1. In particular, the equality holds ex- 
actly in the aforementioned case of 7a = Pa = 0, thanks 
to some nontrivial simplifications of the expressions when 
T[n] is thermal. 



C. Two-time correlators 

We now turn to the problem of the steady state optical 
emission spectrum, that consists in computing two-time 
correlators of the type {c'{Q)c{t)) with c = a, a. We can 
link the two time correlators to the quantities derived in 
the previous sections following an implementation of the 
quantum regression theorem that relies explicitly on the 
density matrix p: 



(ct(0)c(r))=^pf,,](r) {l\c\k) 



(49) 



k,l 



where = (c'^dO {k\){T)) is in the Schrodinger pic- 

ture, where the states evolve and operators have their 
steady state values [S3]. The indices k, I go through all 
the states in the system Hilbcrt space {k — (fci, ^2) with 
ki = 0,l,...,nmax for the photons and ^2 — 0,1 for 
the emitter). The elements Pjj,.;] follow the same master 
equation as the density matrix elements P[k-i]- Since sim- 
ilar approximations can also be naturally implemented, 
this will allow us to provide closed-form solutions for the 
two-time correlators, as is detailed in Appendix |D] We 
give here the main lines of the derivations and introduce 
the key quantities that lead to the final result. We sin- 
gle out, again, only the nonzero elements. For two-time 
correlators, they can be gathered in four functions of n: 

= p[„,»;„-i,»] , i = 0, 1, n>0, (50a) 
<3N = P[„,i;".o] ' ">0' (50b) 

= P[„,0;n-2,1] - n>l. (50c) 

The two-times correlators follow from these quantities 
(with c = a, a respectively) as: 



(at(0)a(r))= VQM 



n=0 



(51a) 



(a^(O)a(T)) = y (V^mS'o[n + 1] + V^Si[n]) . (51b) 



n=0 



Each term n of these sums accounts for the transitions 
between the rungs n + 1 and n, as in the case of sponta- 
neous emission |16j . the first rung, being given by n = 0. 
The equations of motion for the quantities in Eqs. ( 50 ) 



are extracted from the master equation (cf. Eq. D5), and 



into a slow photonic dynamics embedded in a superop- 
erator I?phot on the one hand, and a fast emitter and 
coupling dynamics on the other hand: 

drSo[n + 1] =I?phot{^o[n + 1]} (52a) 

+ ig{VnV[n + 1] - ^/n + lQ[n]) , 
drSiln] =I?phot{^iW} (52b) 

-I- PaSo[n] - JaSl[n] 

- ig{Vn + lV[n + 1] - y/TiQ[n]) , 
drQ[n] =I?phot{QN} (52c) 



+ ig{y/nSi[n] — \/n + lSo[n + 1]) , 
drV[Ti + 1] =Vphot{V[n + 1]} (52d) 

+ ig{y/nSo[n + 1] - Vn+ lSi[n]) , 

for n > 1. After some long, but straightforward alge- 
bra, we can express S'o,i['^] and Q[n] in terms of po,iH 
and qi[n], which, in turn, are expressed in terms of the 
statistics T[n\. This can be done for arbitrary parame- 
ters, including detuning. However, simple expressions are 
possible only at resonance, where we can write Eq. (49) 
as: 



(c^(O)c(r)) =E^ + n, 



00 

■El 

n=0 



C^[n]e-''^°["l^e-(3r"+^^)^/'' + R.s.i.] . (53) 

"R.s.i" stands for "Rabi sign inversion" and is the oper- 
ation that consists in changing the sign of i?i,o keeping 
all other quantities the same. The first term E'^ that 
factors out of the sum, is independent of r, due to the 
approximation of very large photonic lifetime. It is dis- 
cussed separately in section [IV C[ The most fundamental 
quantities that arise in the above treatment are the nth 
manifold inner and outer (half) Rabi frequencies: 

RojM = ^g'{V^^±Vnr - i^^^^Y- (54) 

The (half) Rabi frequency of the linear regime, i?o is 
recovered as the particular case n = with 



i?o = i?o[0] =i?i[0] = Jg 



70 



(55) 



The coefficient Cf[n] (id. for CqIu]) are complex quan- 
tities in general that we decompose into their real and 
imaginary part as 



can also be, like for the single time dynamics, separated 



C^[n] = Lf[n] + iK^[n] 



(56) 
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They are a function of Ri[n] (and Co of Ro[n]). 

As stated previously, the coefficients Cfo[n] are writ- 
ten in terms of the system parameters and the steady 
state photon distribution T[n] only. Their general full 
expression are too long to be given. Only for the simplest 
case of 7cr, 70 = and c — a, the expressions simplify 
sufficiently to be reproduced here: 



— T[nJ H 



-Tin 



where: 



(f)%5'(l + n) 



PS 



+ 



1 + n) 4i?i,o[n] 



(57) 



(58a) 



(f ) -g^(l + nT2Vn(l + n)) 
p2 + 8.g2(i + „) 



3^^)x 



(58b) 



2ff2(^n(l + n)±n) + P^jy^njl + n) T n) 
4(8g2n + P2)(4g4 + 4g2p2(i + 2n) + P^) 



where the notations ± and =p associate / to the upper 
sign and O to the lower one. 



IV. MOLLOW TRIPLET UNDER 
INCOHERENT EXCITATION 

The luminescence spectra can, like all other quanti- 
ties, be obtained "exactly" through numerical compu- 
tations [T^. In our case, computing Eqs. (52) and 
applying the definition of Eq. 



for the power spec- 
tra, one arrives to results such as those shown in Fig. [7j 
where PL spectra are computed both for the cavity and 
the emitter, as a function of pumping. In insets (a-f), 
we select various snapshots plotted in log-scale that il- 
lustrate the regimes discussed previously, namely: (a) 
shows a case from the linear regime, where only the first 
rung of the Jaynes-Cummings ladder is occupied, the 
system being otherwise in vacuum. This yields the vac- 
uum Rabi splitting. The differences in lineshapes in this 
regime have been amply discussed elsewhere (53] . In (b) , 
one is in the quantum nonlinear regime, where the spec- 
trum has a complex structure featuring the many peaks 
that arise from transitions between the lower rungs of the 
Jaynes-Cummings ladder. In (c), one leaves the quantum 
regime, with a collapse of the quantum nonlinear peaks 
that "melt" to form an emerging structure of much re- 
duced complexity, namely, a triplet, as seen in (d) where 
one enters the lasing regime that develops in (e) and is 
fully formed in (f). The triplet is particularly evident 
in the emitter power spectrum, where it is neatly visible 
also in a linear scale. Quenching is not shown explicitly 
in this plot but is elsewhere ^6]. 

A remarkable feature of lasing in strong-coupling arises 
in the form of a scattering peak, well known from Mol- 
low's results (cf. Eq. (13l) where it is due to the driving 
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FIG. 7: (Color online) Cavity (left blue) and emitter (right 
pink) spectra of emission computed numerically. As a func- 
tion of pumping, a transition from (a) the quantum linear 
regime to (f) lasing can be followed passing through (b) 
the quantum nonlinear regime that (c) melts into (d) and 
(e) structures of much reduced complexity, namely, triplets. 
A Mollow triplet is neatly visible in the emitter spectrum, 
whereas the cavity gives rise to a single narrowing Hue. Note 
that in the cuts (a-f), spectra are displayed in log-scale. Pa- 
rameters are 7^ — O.lg, 70- — 0.00334(;, = 7(/> = A = 0. 



laser scattering photons off the atom. This 5 peak also 
forms in the cavity QED version of this problem, as seen 
in Fig. [7](f ) on the emitter spectrum as the very sharp 
peak sitting on top of the incoherent Mollow triplet. Note 
that this peak arises from a numerical computation. It 
is a close counterpart of the Rayleigh scattering peak of 
resonance fluorescence, although this is the coherent field 
grown by the emitter itself in the lasing process that is 
the source of the scattered photons, on the very emitter 
that created them in the first place. More interestingly, 
this peak which, in the lasing regime, is maximum and 
narrower (with the same linewidth as the cavity), forms 
smoothly from a similar depletion when approaching the 
lasing regime from below, as seen in Fig. 7|^d), where an 
equally narrow absorption peak is carved in the emerg- 
ing triplet. As opposed to scattering, in this case, the 
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sj^Ci^ ^ 'V^^^'f"^' 5.9 6.18 6.47 6.77 7.09 7.43 

FIG. 8: (Color online) Lasing as a condensation of dressed states: transitions energies between the dressed states when solving 
numerically the dynamics exactly are shown. The color code (online) has blue shades corresponding to a positive weight of 
the transition in the cavity emission and red shades corresponding to a negative weight. In (a) all regimes are shown over all 
energies, in (b) a close-up of (a) is given for frequencies \ujp\ < g lying between the vacuum Rabi splitting, in (c) a close-up of (b) 
is given in the lasing regime, showing an extremely complex structure of the exact solution obtained in the full-quantization 
picture, although the final result washes out completely most of this underlying information, to provide only the single narrow 
line of a lasing system. 



cavity is coherently and resonantly absorbing excitations 
from the emitter (this depletion is pinned at the cavity 
energy). As the field is building its coherence, it sucks en- 
ergy very efficiently from its source, until a point, shown 
in Fig. [rjje), where the cavity does not require such a 
coherent absorption to keep building in intensity. This 
marks the transition between the point where the system 
is building its coherent field to the one where it is fully 
formed and acting back on the emitter. 

These results are certainly beautiful (an animation 
of this transition is available in the supplementary ma- 
terial of Ref. |55| ) but obtaining them numerically is 
an intensive computer task since the Hilbert space be- 
comes very large, whereas the final output certainly looks 
amenable to a simple analytical description. Solving the 
system numerically nevertheless gives access to every sin- 



gle transition that occur in the many rungs of the Jaynes- 
Cummings ladder through the dressed state decomposi- 
tion of Eq. ([s]) . This yields a surprisingly complex struc- 
ture, shown in Fig. [Sj where we plot the positions ujp 
where the system emits, weighted by the intensity Lp 
of the transition such that resonances disappear with 
vanishing intensities. Fig. |8|^a) gives an overall picture 
while a zoom of the transitions lying between the Rabi 
doublet is given in Fig. [Sjjb), showing an emerging and 
intricate structure, further zoomed in (c) in the lasing 
region. The inner peaks form "bubbles" that open and 
collapse around the origin, where lies the lasing mode. 
Such behaviours of the dressed states also appear in sim- 
pler systems such as two coupled two-level emitters in- 
coherently pumped, which can be solved fully analyti- 
cally [56^. The bubbles formation result from a complex 
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interplay between pumping and decay, which open new 
channels of coherence flow in the system. Figure [sjja), 
shows clearly the satellite peaks of the Mollow triplet. 
Although the lines are neatly splitted the one from the 
other, their increasing broadening allows the formation 
of a smooth spectral shape in the lasing regime, that one 
can follows with the naked eye from the "melting" of the 
quantized structure, as shown in Fig. [t] Fig. |8ja) also 
shows how the inner peaks ultimately all converge at the 
origin, thereby forming the lasing mode. In the lasing 
in strong coupling scenario, lasing can thus be seen as a 
Bose condensate of the dressed states [57H55] . It is fasci- 
nating to follow the formation of a coherent and classical 
field from a fully quantized picture, but this brings little 
insights into the actual phenomenon. Beside hinting at 
its underlying richness and complexity. Fig. [8] essentially 
shows that a complete and fully quantized description of 
a system that is behaving basically classically is hope- 
lessly complicated, keeping track of a huge amount of 
irrelevant details, while the behaviour is well accounted 
for by a few macroscopic degrees of freedom, such as an 
intensity and a off-diagonal coherence element (a). 
Fig. [7] thus illustrates, in one of the most fundamental 
system of quantum optics, the breakdown of the quan- 
tum picture in a quantum to classical transition. Even 
in the simplest and exactly solvable system, it is diffi- 
cult to read much, and we surmise that the condensation 
of dressed states in the lasing process is out of reach of 
the present understanding of dissipative quantum optics, 
calling for a framework such as that developed for con- 
servative systems [5DH5^ . 

Since the intricate patterns of Fig. |8] occur at a dif- 
ferent energy scale than that of the observables and do 
not show up in the spectra, one can hope in the wake 
of the excellent approximations derived previously to get 
similar analytical results also in the spectral domain. In 
the following, we derive such an approximate description 
of the exact picture presented above |21J, allowing us to 
read the essential physics of this transition in the lasing 
regime. 



A. Spectral decomposition 

The Fourier transform in Eq. (|6| of the two-time corre- 
lators ( 53 ) provides an approximated expression for the 



where we have introduced the positions and broadenings: 



spectrum of emission 



TT ^ — ' 



Lf[nr-M_Kf[n]iu;-u:,[n]) , 



R.s.i. 



R.s.i. 



wi.oN 



5R(i?LoW). 
SFo- + 70 



23(i?i,oN), 



(60a) 
(60b) 



so that the optical spectrum is composed of a series of 
Lorentzian lines at frequencies ujp with linewidths 7p and 
weight Lp plus interference terms weighted by Kp. These 
lines arise from transitions between rungs n + 1 and n of 
the Jaynes-Cummings ladder. The spectra are normal- 
ized to unity, therefore, the incoherent part of the spectra 
(second and third lines) is normalized to 1 — di{E'')/nc- 
Each transition can exhibit weak or strong coupling, that 
is, the rungs being split or not into dressed states, sim- 
ilarly to the case with no pumping 16j. When split, all 
peaks broadenings are the same, regardless of the rung: 

71 W = 70 W = (sr,^ + 7(/')/2- 

In Fi g. [9} we see two examples of spectra computed as 
in Eq. (|59[) with T[n] taken as Poissonian. The first case, 
with 7a = 0.1(7, corresponds to a very good cavity well in 
the strong coupling regime with the emitter, such as is 
realized in circuit QED. The second case, with 7^ = O.Sg, 
corresponds to a less favourable situation representative 
of the state of the art of systems such as quantum dots 
in microcavities ^HJ". In both plots, the pump increases 
from the top to the bottom curves and features the quan- 
tum to classical crossover. The transitions between the 
Jaynes-Cummings rungs are first resolved individually, at 
low pumping, and then merge into a Mollow triplet. The 
approximated spectrum of emission differs from the nu- 
merical exact result at low pumpings, since the assump- 
tion of much larger pumping than decay does not apply 
here. The regime of validity for our approximated spec- 
trum and that for the observation of the Mollow triplet 
with incoherent excitation is: 

7a , 7o- , 7* < 5 < -Po- < • (61) 

Note that out of resonance, one must consider (7eff[l] in- 
stead of g in Eq. (61 1 . As the position of the peaks is still 



well approximated, however, this approximation provides 
an instructive and physically transparent picture of the 
Mollow triplet formation. We analyse these peak posi- 
tions in more details in the next subsection. 



B. Peak positions 

Let us recall the expression of the Jaynes-Cummings 
quadruplets positions in the spontaneous regime, in order 
to appreciate better the features brought by the incoher- 
ent pump [TB] : 



la - la 



(59) 



± \ g'^n- 



la - la 



(62) 
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FIG. 9: (Color online) Approximated spectra of emission, 
Eq. ( |59| l, as pumping is increased and the system undergoes 
the quantum-to-classical transition. The plots correspond to 
two different cavity lifetimes {■ja ~ O.lg and 0.5g). Parame- 
ters, not specified in the plots, are Pa = 7o- = 7</, = A = 0. 
In inset, the comparison between the approximated-analytical 
(thin blue) and exact-numerical (thick black) spectra for the 
lowest pump Pa = O.lg. Although the approximation breaks 
in the quantum regime, it shows how the MoUow triplet is 
formed. 



These branches are plotted as thin dashed lines in Fig. [TO] 
(only the positive frequency ones) as a function of their 
decoherence rate T = 7a — 7cr. They can be exactly 
mapped with transitions between rungs in the dissipative 
Jaynes-Cummings ladder, exhibiting three regimes de- 
pending on their splittings: (a) at low decoherence, both 
rungs — initial and final — are split (in the strong coupling 
regime) and we have four peaks (two shown in the fig- 
ure), (b) When the lower rung enters weak coupling and 
recovers the bare modes, only two transitions (one in the 
figure) can be seen from the upper, still split, rung, (c) 
At large decoherence, none of the rungs are split and all 
transitions are at the origin. 

In the same Fig. [TOj we have plotted with thick lines 
the positions in the presence of pumping, Eq. (60a), 
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FIG. 10: Transition energies (positions of the spectral peaks) 
as a function of decoherence |r|/(4(;) from rungs n = 1, 2, 3, 4 
(from dark to light lines) for two cases: spontaneous emission 
(r = 7a — 7ct, Eq. (621, thin dashed lines) and in the pres- 
ence of pumping (T — — 70, Eq. ( |60a[ | thick lines). The 
line starting a,t uj — g corresponds to the linear regime or 
first-rung-to-vacuum transition (5R(-Ro), thick black line) and 
is common to both cases. The curves starting above (resp. 
below) 1 correspond to the outer (inner) transitions. 



as a function of the corresponding decoherence rate, 
F = Per + ja- The positions, paired around the centre 
(still only the positive half is shown), also depend on 
the rung. In the linear regime, n = 0, these are simply 
the two Rabi transitions ±3?(i?o), with the same struc- 
ture as two coupled harmonic oscillators or the sponta- 
neous regime that we have just described. However, for 
n > outer and inner peaks have an interesting non- 
trivial behaviour. The inner peaks close at low F, before 
the Rabi transitions do, forming, in the lasing or clas- 
sical regime (ria ^ 1), the central peak of the MoUow 
triplet. The outer peaks close much after the correspond- 
ing Rabi peaks do, grouping to form the two side bands 



(r./4)' 



of the MoUow triplet with ~ 3fi^/4ng2 

(see Fig.|9]). 

In the presence of pumping, it is no longer possible 
to map exactly the position of the peaks composing the 
spectra with transitions from the Jaynes-Cummings lad- 
der. The effect of the incoherent pump, before being so 
strong to close the dressed state splitting of the Jaynes- 
Cummings rungs, is to make it homogeneous throughout 
the ladder. Then, the inner peaks do not close because 
the rungs enter the weak coupling regime but rather be- 
cause the transitions coincide with the cavity frequency. 
In any case, the Jaynes-Cummings structure is very much 
distorted and the transitions between its rungs very much 
mixed, in such a way that it is no longer possible to re- 
construct the ladder. Furthermore, in our derivation of 



Eq. ( 54 1 we have neglected the interplay between the cav- 



ity and emitter dynamics (by assuming them to have very 
different time scales) and, therefore, Fig.[9]shows the iso- 
lated effect of the incoherent pump in contrast with the 
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complex features we observed in Fig. [8] 



C. Elastic scattering E'^ 

It is possible to give the general expression (with de- 
tuning, dephasing and decay) for the elastic scattering 
contribution to the spectra, for both emitter and cavity 
emission: 



4.9 



■70 



E 



(63a) 



7^(Ka + 2r^)Vl + nqi 


n 


+ 1]^ 


- Paina- 2T„)y/nqi[n] 




-i 


-4r2 


+ 4r,K,(2n + l) 



K(i;°) = y . 

fo [l + (p^)2j + ATI + 4r.AJa(2n + 1) 
X {7a(Ka(4n + 1) + 2r,)(l + n)T[n + 1] 
+ P^(K„(4n + 3) + 2r„)nT[n]} . (63b) 

Note that U{E°-) > in all cases, while ^{E") can be 
negative, as we pointed out in the discussion of Fig. [7] 
The expressions (63) support our qualitative discussions. 



In the cavity emission, the scattering peak is simply the 
direct emission of photons through the cavity mode and 
therefore, is always positive. For the emitter, however, 
these coherent cavity photons must "convert" into mate- 
rial excitations before being emitted, which implies the 
possibility of interferences that can be positive or nega- 
tive. This change of sign is even more clear when 7o- = 
where E"^ < if Po- < V^ff- At low pumping, this carves 
a "hole" in the spectra pinned at the cavity frequency. 
For Pa > \/2g, the delta peak is positive. In the case 
of the cavity, most of the emission comes from it, and 
not the de-excitation of the dressed states, that is, from 
the photons that undergo an efficient interaction with the 
emitter. 

When the two-level system is not sufficiently popu- 
lated/inverted (at low Pg.), the cavity has a smaller elas- 
tic scattering contribution. The emitter sees an inter- 



ference hole being carved in its spectrum, reminiscent 
of a Fano resonance. The emitter represents an efficient 
pumping for the cavity in a linear way. As we already 
mentioned, the cavity is sucking the cavity photons (at 
uja) out of the emitter. Although at very low pumping 
(Pg. « 7^,) these approximations are not valid (a delta 
function weighted negatively means a negative spectrum) 
and do not provide a quantitative agreement with the 
numerical results, they are interesting to understand the 
qualitative features that are obtained numerically. 

When the two-level system is saturated in the self- 
quenching regime (at high Pa-), it is mainly in its excited 
state. The cavity spectrum tends towards the bare cav- 
ity emission (thermal regime). In the exact solution, this 
is a Lorentzian with the cavity linewidth 7^. In our ap- 
proximated solution, the cavity bare emission is a delta 
function and the two-level system emission loses com- 
pletely the elastic scattering component because there is 
weak interaction with the cavity which provides it. 

In the lasing region, 1 < P^r < 20g, the triplet appears 
in both channels of emission: cavity and emitter. The 
emitter grows some small positive elastic scattering com- 
ponent on top of the triplet, as a collateral effect from 
the strong interaction with the cavity (proportional to 

qiN)- 



D. Semiclassical approximation 

In the lasing regime, where the n and n + 1 rungs that 
are close to (having rij, ^ 1) have similar splittings 
and the cavity field becomes Poissonian (coherent), sim- 
ple expressions — in regard to the complexity of the un- 
derlying machinery to obtain them — arise for the spectra 
of emission [21^ . The variable n becomes continuous as 
compared to the large intensity tIq and, given that the 
distribution is Poissonian (peaked at the mean value), 
we can consider only the case n = in Eq. ( 59 ) . The 
integral over the distribution T[n] simplifies to 1. In this 
regime, the inner peaks have collapsed into the centre 



while the outer remain splitted. Substituting Ua (and 
for the normalization) from Eq. ( 59 ) , we obtain the final 
expression for the emitter spectrum: 



^2PaTJSTl 



2Pa 



_^ 2 

27r / r„+7^ N2 



- 70 + '*.)[(-2P. + TaYnl + ((3r, + 7^)2 + 4(r, - 2P<,)ac,)c^2 + 4^4] { 4P>,(3r, +7^-1- n^) 
f4r<,(7^ + 2K^) + (70 + K<T)(70 + 3K^)]+4P<,K^w2-(r^+7^ + K„)[r2(3r^-H70-H2K<,) + (r^-70)w2]| 



(64) 
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FIG. 11: (Color online) Comparison between the position of 
the side peaks (5R(_Ro), in solid blue) and those observed in 
the emitter spectrum (cjobs, in solid purple, above the filling), 
as a function of pumping. In dashed is shown the neck of the 
side peak so that the filled area delimits approximately its 
half-width. 



Similarly to the coherent excitation case, it is composed 
of an elastic scattering term (delta peak) , a central peak 
(a Lorentzian peak with FWHM To- +7,/)) and two side 
bands. When splitted, these have a FWHM (SF^ +70)/2 
and positions given by ojoi'T-a] ~ 5J(i?o) where 



This is the Mollow splitting in the case of incoherent 
excitation, analogous to i?L in Eq. (11). It is plotted 
in Fig. 
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Contrary to the laser excitation, this split- 
ting can now close due to the pumping intensity P„. In 
Fig. [TT} we show the domains where the Mollow triplet is 
clearly resolved. This naturally requires that the system 
is able to enter the regime of lasing in strong-coupling, 
which starts at figures of about 7a/5 ~ 0.1. This is the 
case shown in the Figure, where we compare the Mol- 
low splitting, Eq. (11), with the observed splitting, rep- 



resented as the shaded area which is delimited by the 
maximum (upper boundary, solid) and the neck (lower 
boundary, dashed) of a side peak. When the Mollow 
splitting is maximum, decoherence has however broad- 
ened so much the satellite peaks that no triplet is ob- 
servable anymore. 

Applying the same procedure, for the cavity emission 
we have Sa{y>) — (^(w), that is, a purely elastic spectrum 
with negligible linewidth. A more accurate approxima- 
tion to the FWHM of the elastic peaks in this regime, 
which reproduces the typical lasing cavity line narrow- 
ing, is given by the broadening 7l: 



7L = -Ig^lalPl 

as derived by Poddubny et al. |17j . 



(66) 
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FIG. 12: (Colour online) Mollow triplet in the direct PL 
spectra from the emitter, Scr{ijj) for ja = O.lg, Pa — 0, 
= 0.00334 and = 7g (ua « 11.6 and n„ « 0.53), for 
various configurations of detuning and dephasing. The exact 
numerical results are plotted with a dark-blue dashed line, 
overlapping almost exactly our analytical formula, in light- 
blue solid. The coherent scattering peak is featured in the 
analytical solution in (a) while only the incoherent spectrum 
is shown in (b) and (c). The dotted line in (a) represents 
the Mollow triplet obtained under coherent excitation, with 
a laser intensity equivalent to the cavity occupation in the 
incoherent case, that is, for JIl = \/n^g and also an equiv- 
alent emitter broadening, 70- + Pa- The two types of 
Mollow lineshapes are clearly different even though their peak 
positions and broadening are equal. Also, as compared to the 
coherent excitation case, the Mollow triplet under incoher- 
ent pumping is a resonant structure that becomes strongly 
asymmetric with detuning. 



In Fig. [T2I we show the excellent agreement between 
the analytical approximation (in dashed lines) and the 
exact numerical computation (solid lines). In Fig. 12 'a), 
the case of resonance — the one of most interest — is where 
the Mollow triplet is best observed. Its analytical ex- 
pression is given by Eq. ( 64 1 , where we also include the 



scattering peak, (both as a result of the numerical pro- 
cedure and approximated by Eq. (66)) seen as a very 



sharp line, which we have truncated in the plot, as it ex- 
tends more than one hundred times higher than is shown. 
In Fig. [l2{^a), we also provide further evidence that the 
Mollow triplet formed under incoherent pumping is of 
a different nature than that formed under coherent ex- 
citation |21| . by superimposing the coherent excitation 
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(a) incoherent excitation (b) coherent excitation 




0.25 



''o 10 20 30 40 10 20 30 40 

PJg Pa/g 

FIG. 13: Contributions to the MoUow triplet under (a) in- 
coherent and (b) coherent excitation for the three types of 
peaks: the central Lorentzian peak (constant at 1/2, solid 
line), the side bands (dashed) and the elastic scattering peak 
(dotted). The parameters are chosen so that the two MoUow 
triplets can be compared on equal grounds, as explained in 
the text, for 7a = O.lg, Pa ~ 0, ja = J4, ~ and A = 0. 



Mollow triplet (dotted line). We take 17l — \/n^g and 
7(T — > 7(T + Per to compare both expressions on equal 
grounds (both lineshapes remain dissimilar even when 
parameters are left completely free). This substitution 
makes both types of triplet share the same position and 
broadening for their three peaks. However their relative 
weight is different, as shown in Fig. [13) These strong 
qualitative departures result in the striking differences 
between the final spectra, although the peaks have iden- 
tical characteristic if taken in isolation. 

When the system is not at resonance, in sharp con- 
trast with the conventional Mollow triplet that retained 
its qualitative features (cf. Fig. [T]), the Mollow triplet 
under incoherent pumping becomes strongly asymmet- 
ric, as it recovers the scenario of an anticrossing of two 
modes. The out-of resonance case is studied in details 
in Fig. [M] Whereas both detuning and dephasing were 
needed to break the symmetry of the conventional Mol- 
low triplet, the one formed under incoherent pumping is 
lost by detuning alone, pumping playing already the role 
of dephasing. Dephasing has otherwise the expected ef- 
fect of smearing out and broadening the spectral features. 
Analytical results can also be given for the non-resonant 
case when A 7^ 0, but, as for the conventional Mollow, 
they are too long to be reasonably written down. We plot 
it on top of the numerical solution in Figs. [T2j[b) and (c) 
where one can see the semiclassical approximation is ex- 
cellent there as well. 



V. SUMMARY AND OUTLOOKS 

We investigated the steady states of the Jaynes- 
Cummings Hamiltonian established under the interplay 
of decay and incoherent pumping (of the emitter), in- 
cluding pure dephasing and detuning for wider generality. 
This is the simplest and most fundamental realization of 
a fully quantized system, realized with atoms, quantum 
dots or superconducting qubits in a cavity. We identify 
five regimes through semi- analytical and approximated 
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FIG. 14: Counterpart of Fig. [T] for the Mollow triplet under 
incoherent excitation. In contrast to the coherent excitation 
case, detuning alone breaks the Mollow triplet. Parameters 
are 7„/(? = 0.1, 7<,/(/ = 0.01, Pa = and P^jg = 7. 



solutions, all confronted to exact numerical solutions. 
These are: i the linear quantum regime, ii the nonlin- 
ear quantum regime, in the nonlinear classical regime, 
or lasing regime, iv the self-quenching regime and v the 
thermal regime. We provided closed-form analytical ex- 
pressions that account for most of theses regimes and a 
simple numerical procedure (solving a set of two coupled 
nonlinear equations) that afford an excellent description 
over the entire range of excitation and all the five regimes 
that we have outlined. This also allows a transparent 
reading of the physics involved, namely, the first regime 
involves the lowest rung of the Jaynes-Cummings ladder 
only and corresponds to spontaneous emission. A linear- 
model (two coupled oscillators) and a truncated Jaynes- 
Cummings model offer two complementary views of this 
regime. The quantum regime is the one where the system 
starts to climb the ladder, requiring a full record of all the 
quantum correlators involved. This is therefore the most 
complicated regime from the point of view of the amount 
of information required to describe it, since no good ap- 
proximation can synthesize the dynamics of a few quanta. 
In very good systems, this manifests spectrally in a com- 
plex structure of peaks at anharmonic frequencies. As 
pumping is further increased, lasing ensues which brings 
back the system to a simple level of description in a semi- 
classical approximation. A single narrowing line in the 
cavity or a variation of the Mollow triplet for the emitter 
describe accurately the system. In the fully quantized 
picture, lasing appears as a condensate of dressed states, 
with a complex pattern that however does not manifest 
in the observable, showing a breakdown of the quantized 
picture in favour of a classical description. 

The Jaynes-Cummings model that arose as a challenge 
for full-field quantization [M] remains to this day a pro- 
ficient source of theoretical and experimental investiga- 
tions into the quantum realm. The realization in the 
laboratory of the nontrivial quantum physics that it cov- 
ers will shed light on foremost issues such as the quan- 
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turn to classical crossover, emergence of coherence, las- 
ing and quantum nonlinearities. A solid understanding 
of the various regimes it realizes may also lead to useful 
devices and applications, from single-photon sources to 
low-threshold lasers. 
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Appendix A: Quantum regression formula for 
coherent excitation 



Two-time correlators of the type (cr^(O)cr(T)) can 
be computed by means of the quantum regression 
formula 16 . Once we find the set of opera- 
tors „} = (J^ (t" (with m, n e {0, 1}) and the 
regression matrix M mn that satisfy Tr[C/m,n}'CO] = 

Sim'nn ^"}",Tr[C/m',n'}0] for a general operator O, 
then the equations of motion for the two-time correla- 
tors (r > 0) in the steady state {t — 0) read: 



The solution is 



E 



9,(O(0)C{„,„}(r)) = M"- (O(0)q™',„'}(T)), 



(Al) 

The corresponding regression matrix M is given, in the 
case of a coherent and classical excitation of the emitter, 
as explained in section In) by: 



Mmn — iAr(m — n) — —(m + n) — — (m - 

run ' 2 ^ '2^ 



Al mn 
1 — m.n 

Al mn 
m , 1 — n 



iQi^lm + 2n{l — m)] , 
-~ifli^[n + 2m(l — n)] 



(A2a) 
(A2b) 

(A2c) 



and zero everywhere else. We concentrate on comput- 
ing g[^\ which corresponds to setting O = and hav- 
ing {m, n} = {0, 1} in Eq. ( Al ). We obtain the following 
differential equation: 



arVL(r) = -MiVi(r)+Ai(at) 



where 



VL(r) = 



(^t(o)a(r)) 
(at(O)at(r)) 

((7t(0)<TV(T)), 



At = iil 




(A3) 



(A4) 



and 



la +14. 

i ' 2 










(A5) 



vl(t) = e 



[vi(0) - M-iAi(at)] + M-iAi(at) , 

(A6) 

in terms of the steady steady valu es {cr ^) and vl(0) = 
[((tV), 0, 0] (initial condition of Eq. (A4)). We also need, 
therefore, to compute the steady state of the system, 




(A7) 



which can be done, again, by means of the general for- 
mula in Eq. (Al). This time we take O = 1 and 
{m, n} = {0, 1}, {1, 0}, {1, 1} and find the equation: 



dtVLL = = Mlvll + Al . 



(A8) 



The solution is — M^^A^. It allows us to simplify 



Eq. ( A6 1 further as 



vl(t) 



-Mlt 



[vi(0)-Ui(at)]+Ui(at), (A9) 



The final explicit solutions for the mean values and cor- 
relators of interest are presented in the main text. 



Appendix B: Derivation of the field correlators and 
density matrix 



The equations of motion of the correlators in Eq. (17) 



can be derived from the master equation ( 16 ) by simply 
applying the general formula (O) = Tr[Op) or from the 
rules given by the quantum regression formula jl6j : 

dtNa[n] = -nTaNa[n] + n^PaNa[n - 1] (Bla) 
+ 2gnN'^M, 

dtN,[n]^-[ra{n-l)+r,]N,[n] (Bib) 
+ Pain^ l)^N,[n - 1] + P,Na[n - 1] - 2gNlM , 

K.M (Bic) 



+ Pan{n ~ l)NUn - 1] + AN:M 
+ g{2N,[n + 1] + nN,[n] - 7V,[n]) , 



dtKM 



f(2n-l) + ^2 



KAn] (Bid) 



■Pan{n-l)N:,[n-l]~ANlM, 



with n > 1 and also Na[0] — 1. In the steady state, 
with Pa = 0, one can further simplify these equations 
and write Naa[n] and N„[n] in terms of A'a[n] to obtain 
Eqs. (IT9|. 



The master equation (16) can also be rewritten in 
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terms of the density matrix elements in Eq. (45) as: 

(B2a) 



dtPo[n] = -[(7a + Pa)n + Pa\po[n] 
+ 7a(" + l)Po["- + 1] +Panpo[n~ 1] 
- PcrPoN + To-PiW - 2.gV"qiN , 

<9tPlN = ~[{-fa + Pa)n + Pa]pi[n] 

+ la{n + l)pi[n + 1] + Panpi[n - 1] 



(B2b) 



dtqi [n] 



7crPiN + PaPQ[n] + 2g^/n + lqi[n + 1] , 

+ 7aV("- + l)»^qi[" + 1] + Pa\/{n- l)nqi[n - 1" 
+70 



(B2c) 



-qiW + ^^"(PoW - Pi["^ - 1]) - Aqi.[n] , 



9tqrW 



(7a + Pa)n ■ 



Fa 

2 



qrN 



(B2d) 



+ 7a \/[n + l)riqr [n + 1] + Pj, ^(n - l)nqi. [n - 1] 
- ^-^^Y^qrW + Aqi[n] . 

At resonance, the real part of the coherence distribution, 
Qr, gets decouple d and v anish es in the steady state. As a 
result, only Eqs. (B2b)-(B2d) need to be solved. When g 



vanishes, qj does not couple the two modes anymore, and 
their statistics become thermal, like in the boson case. 
Through the off-diagonal elements qi, the photon density 
matrix can vary between Poissonian, thermal (superpois- 
sonian) and subpoissonian distributions. 

In order to solve these equations in the steady state, we 
neglect the photonic dynamics (as explained in the main 
text) and further substitute qr [n] — [2A/(rcr+70)]qi[fT.] in 
the equation for q;, and write everywhere po[?^] as T[n] — 
Pi[n] and pi [n + 1] as T[n + 1] — po[n + 1], so that qr[n-], 
Po[n], pi[n+ 1] do not appear explicitly in the remaining 
three equations. Then, the equations read in matricial 
form: 



9tUo[n] = -Mo[n]uo[n] + Ao[n] 



(B3) 



with: 



Uo 71 



Mo[7i] = 



Po["^ 






Pi 




, Ao[ri 


qi {n 













( 









2gVn+ 1 



-g^n + 1 g^n+1 



-2gy/n + 1 



(B4) 



The solution in the steady state is uo[n] = 
(Mq [n]) ~"'"Ao[n] which gives the result of the text, 
Eq. (47 1. Note that this solution is only exact in the 
case 7„ = Pa = 0. 



Appendix C: Perturbative regime of interaction in 
the limit of wealt coupling 



Eq. ( 20 1 can be solved exactly as a series Taylor expan- 



sion on the pumping Pc. For this, we rewrite Eq. (20) 
in terms of the fraction F[n\ ~ Na\n + IJ/A^aM = 
"a5^"+^V5^"^ (for " > 1): 



F[n-l]{F[n] 



where 

■^n = 

B 

a, 



27a 



+ n7a 
1 



'17a 



2P„ 



nP„ 



r^r + (n - l)7o 



(CI) 

(C2a) 
1, (C2b) 
(C2c) 



All quantities can be expanded, or assumed to have a 
solution in the case of P[n], in power series of Py. 

„ oo „ oo 

fc=0 " fc=0 



fc=0 " 
oo 

PM = 5^A.[n]P;. 



(C3b) 



A;=0 



A key feature of this expansion is that /3o['^] = and 
/o[n] = for all n (we recall the linear behaviour of Ua 
at low pumping). With this considerations, Eq. (CI ) now 
reads 



oo oo 



l](/,[n] +afc[n])Pr^ = Y,mPl ■ 

q=l k=l t=l 

(C4) 

We further change the sum index q, on the left hand side 
of the equation, for t = q + so that we can get rid of 
it and of the pump dependence: 

t 

E ft-k[n - l]{fk[n] + ak[n])P^+'' = I3t[n] . (C5) 

k=l 

The exact solution can be found exactly and recurrently 



fi[n] 
ft[n] 



I3i[n + 1] 
ao[n + 1] ' 

/?2[n + l] -/iM(/i[" + l]+Qi[»+l]) 
ao[n + 1] 



(C6a) 



(C6b) 



f3t[n + 1] - ft-i[n] EtJiifkln + 1] + ak[n + 1]) 



aQ[n - 



(C6c) 
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This method is only useful in practical terms when the 
effect of the coupling is perturbative, that is, at very low 
pumping or in the weak coupling regime, where only a 
few terms of the expansion are needed. Otherwise, in 
order to reproduce non perturbative effects such as the 
transition into lasing, the sum should be performed to all 
orders of P„, which might not be practical numerically. 



Appendix D: Derivation of the two-time correlators 

We can obtain the elements p%.i{t) needed to compute 
the two-time correlators in Eq. ( [49| ) , in an equivalent way 
P[k;i]j &s they follow the same master equation |53j : 
9tP[k;i] = ^ {h-i] P[i;i]- That is, we can solve 



SrPffe;,] (i + r) - ^ X r ^ . pf,^^.] (t + r) . (Dl) 

The initial conditions are given in terms of the steady 
state density matrix: 

Pffc;Z](0) = {c\\l) {km) =Y.Pik-M («l 10 . (D2) 



Let us be more specific by writing the formulas for the 
two correlators of interest. For the emitter spectra, we 
have (^1, /2I cr |fci, A:2) = ^;i,fci(5;2,o^fc2,i' which gives: 



(D3) 



n=0 



For the cavity spectra, we have (^i, Z21 1 l^i, ^2) 
\/ki5i^^k^5i^M-i^ which gives: 



ri=0 'i=0.1 



p'^ is obtained by solving the master Eq. (Dl| in both 



cases, but the initial conditions that follow in each case 



from Eq. (D2), are different. For the emitter correlator, 
they are pf„_i.„j](0) = P[n,i-mS]5j,o, while for the photon, 

/'fn,«;mj](0) = V?" + lp[n,i:m+i j] ■ 1^ the Same way as 
when solving the steady state distributions, we write the 
master equation for p? • only for the elements that 

[71, i, 772.^ J 

will be different from zero during the evolution with r. 
We have to include all the elements that are nonzero in 
the initial condition plus those that are linked to them. 
One can check that the nonzero elements are the same for 
the initial conditions of both the cavity and the emitter 
correlators, those defined in Eq. (50). They follow the 



equations: 



ria+Pa 



(2n-l) + PJ5oN 



(D5a) 



+ 7av/n(n + l)S'o[n + 1] + Pa\/n(n - l)S'o[n - 1] 
- PaSa[n] + -laSM + ig{V^{^V[n] - y/nQ[n - 1]) 

na+Pa 



drSi[n] 



-{2n-l) + Pa]Si[n] 



(D5b) 



+ -ia^/n{n + l)Si [n + 1] + Pa^/n{n-l)Si [n - 1] 
+ PaSa[n] - laSi[n] - ig{^/n + lV[n + 1] - ^/^Q[n]) 



drQ[n] = - [(7a + Pa)n + Pa\ Q[n] 
+ la{n + l)Q[n + 1] + PanQ[n - 1] 
'Ta +70 1 



(D5c) 



(D5d) 



+ ig{\/nSi[n] - ^/n + lSa[n + 1]) , 
drV[n\ = - [(7, + Pa)n - 7j V[n] 
+ -^a^/{n-l){n + l)V[n + 1] + Pa^/n{n ^ 2)V[n - 1] 



+ ig{\/n - 15*0 M — y/nSi[n - 1]) . 

As explained in the main text, we solve these equations 
by neglecting the very slow photonic dynamics. We then 
define the steady state and slow varying function 



X[n] = S'o[n](0) + 5i[n](0), 



(D6) 



and substitute 5*0 [«] = X[n] — 5*1 [n] and Si[n + 1] — 
1] — S'o[n.+ 1], so that we can rewrite the equations 
in a matricial form (IBSl: 



drni[n]{T) = -Mi[?i]ui[7i](r) + Ai[n] 



(D7) 



with: 



Q[n] 
\V[n+l] 



, Ai[n] 



PaX[n] 





and Mi[n] = 







igy/n + 1 
\ —igy/n ig\Jn + 1 



ig\/ n + 1 —ig^fn \ 

—ig^fn ig\/n+l 

igy/n 



^^-lA 



+ ^A/ 



(D8) 



The solution in the steady state is 



ui [n](r) = e-^i["l^(ui[n](0) - iM,[n])-' A,[n]) 

+ (Mi[n])^iAi[n]. (D9) 



22 



For the emitter correlator, the initial condition derived in Eq. (53), E'^, which turns into a delta function at the 



from Eq. (D2) reads: 



ui[n](0) = 



V 








(DIO) 



cavity frequency in the frequency domain, Eq. ( 59 1 



For the emitter, the intensity of this contribution is 
found from summing the third element of Bi[n], for all 
rungs n>0: E'^ = X]^o(^iN)3- Fo^' the cavity case, it 
is found as E'' = J2n=o [V?H^(Bi[n])i + Vn(Bi[n])2 



which implies substituting X[n\ = r ^+7^, ^" *jQiM 
Ai[rt]. The final correlator is found by taking the third 
element of the vectorial solution and summing contribu- 
tions from all rungs: (cr^ (O)g(T)) = X]5^Lo('^i["])3' which 



corresponds to Eq. (51a) 



For the photonic correlator, the initial condition is: 



UiW(0) = 



+ Ipo [n + 1] 



2A 



(Dll) 



which implies substituting X[n\ — Y^T[n] in Ai[n]. 
The final correlator is found as: (a^(O)a(r)) = 

J2^=o Vn + l{ui[n\)i + ^/n{ui[n\)2 , which corresponds 



to Eq. (51b) 



1. Elastic scattering term 



The second line in Eq. (D9), 



Bi[n] = (Mi[n])-iAi[n] 



(D12) 



is independent of r, due to the approximation of infi- 
nite lifetime 7^ ~ 0. If ja was of the order of g,Tcy, we 
could not have assumed X[n] to be r-independent and, 
therefore, nor Bi[n]. In the regime where our approxi- 
mation is valid, this term gives rise to the first element 



2. First rung 

In principle, we must solve separately the Rabi dynam- 
ics of the first rung with the ground state, n — 0, having 
a 2 X 2 system {Si[0] = and V[0] = 0): 



ui[0] 



[Soil] 



Ai[0] 







Mi[0] 



I'd 



For the emitter, the initial conditions are: 



ui[0](0) = 



In the photonic case, the initial conditions are: 




(D13) 



(D14) 



ui[0](0) 




Po[l] 



(D15) 



One can check that the solution for n = is finally the 
same as taking the limit n — ^ in the expressions we 
obtained for n > 1. Therefore, there is no further need 
of separating the Rabi from the rest of rungs in our ex- 
pressions, although it will find simpler expressions for the 
quantities of interest and we may point them out. For 
instance, if 7^ = 0, in both cases, Ai[0] = and the 
solution simplifies to: Ui[0](t) = e-^i[°l^ui[0](0). 
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